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This paper presents an analytical review of the transient events relating 
to charge, current, voltage, and energy upon disconnection of inductive 
loads in switching circuits. Based on good agreement between theory and 
experiment, the analysis is found applicable to a number of long-standing 
relay and switching circuit problems involving radiation interference, 
voltage breakdowns in various elements of the contact circuit, and contact 
erosion. 

I. INTRODUCTION 

The disconnection of a relay in a switching circuit involves a number 
of time-variant design problems for which a quantitative treatment 
is needed. These mainly concern (i) the release time of the structure, 
or (it) the transient current and voltage surges which are a potential 
cause of oscillations and breakdowns in the circuitry. The subject of 
release time is chiefly related to a change in the core's magnetization 
from its operated to its residual value, i.e., a very substantial change 
in core flux occurring over a comparatively long time, typically several 
milliseconds. The transient circuit surges tend to develop much earlier 
in the interval after disconnection, often in times from about 10" 
to 10~ 4 second during which the core's magnetization has not yet 
had time to change appreciably. Thus, while both phenomena are 
related to the core's demagnetization curve, their performance tends 
to be controlled by different features of this nonlinear characteristic. 
The subject of release time, including the design of slow release relays, 
has been on a firm quantitative basis for many years, 1-4 based on 
relationships involving the whole demagnetization curve. However, 
only qualitative or empirical treatments have so far been available 
to describe the transient surges of current and voltage. An analytical 
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description of these surges is indispensable to the understanding of 
voltage breakdowns or radiation interference in the contact wiring; 
and a renewed interest in these latter subjects has prompted their 
further study as outlined in this report. 

The surges principally result from the tendency for the inductance 
to maintain current flow after the main current path has been opened; 
current is then diverted into the circuit capacitances, temporarily 
charging them to high voltages, even to several hundred times the 
original battery voltage (when the capacitance is very small). Previous 
analyses based on a conventional lumped-element L, R, C treatment 
have not agreed well with experiment, attributed mainly to difficulties 
in suitably defining the load's nonlinear inductance and to the neglect 
of eddy current losses in the core. In attempting to recognize these 
in the present analysis, some simplifications were possible by assuming 
an effective lumped-element value of inductance which corresponds to 
the differential permeance of the load's magnetic circuit in the operated 
region. This concept, and an allowance for eddy current losses, give 
results which agree well with experiment over several decades of vari- 
ation in the basic parameters. Over the three decades of principal 
interest when studying breakdown voltage, this effective inductance 
is found to be essentially constant, being basically determined by the 
differential permeance of the core's demagnetization curve at the 
operated point. Beyond this region, where pulses are slower, the use 
of a progressively increasing permeance, which recognizes the wider 
range of the demagnetization curve that comes into play, helps improve 
the accuracy. The study covered in this report deals mainly with the 
region where the faster pulses predominate. 

1.1 Outline of Analysis Procedure 

Expressions are first obtained for the general case of a relay winding 
having both a short-circuited secondary and an RC network, as com- 
monly used for contact protection, bridging the main winding. By 
assigning suitable values to each of these terms, the results can then 
be simplified to cover a range of practical situations such as the ordinary 
unprotected relay coil, the coil bridged by either capacitor or resistor 
alone, and the air core coil. Of these, the most common case, that of 
a relay coil shunted only by a capacitor (representing distributed 
capacitance, wiring capacitance, etc.), is selected for detailed study. 
Based on good experimental agreement for a wide range of test cases, 
the results are then applied to a number of special relay and switching 
circuit design problems. 
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The surges to be studied are illustrated by the curves in Fig. 1. 
This view shows a typical relay switching circuit in schematic form, 
and the transient current and voltage surges when contact X is opened. 
A major objective is to describe these surges analytically. The initial 
phase of the analysis is to derive expressions for current and voltage, 
as a function of time, when the load is disconnected by an "ideal" 
(i.e., nonbreakdown) contact. The smooth curves in Fig. 1 show typical 
results when a telephone relay is controlled by a (nonarcing) mercury 
contact. They have practical interest whenever the switching is done 
by a vacuum tube, a transistor, or a mercury contact which is well 
enough insulated to avoid breakdowns. Such curves also determine 
the conditions for breakdown when apparatus, wiring, or contacts 
of known dielectric strength are connected in the circuit. For example, 
the irregular curves in Fig. 1 show the multiple breakdowns that occur 
across a typical relay contact in air, as its surfaces separate. This 
result was obtained by replacing the previous mercury contact by an 
ordinary mechanical contact (on a wire-spring relay) and repeating 
the measurement as a double exposure on the same film. Such break- 
downs, sometimes termed "B-type transients" or "showering arcs," 
represent successive charging of local capacitance to the dielectric 
breakdown of the contact's gap, its discharge through the local contact 
circuit, and a renewed charging cycle. They have been described 
qualitatively in considerable detail, 5-8 and there have also been some 
efforts at a quantitative treatment. 9,10 Each successive charging cycle 
is the beginning of a new "ideal" contact surge; hence the definition 
of this curve facilitates the study of such breakdowns. 

A second phase of the analysis examines these surges from the stand- 
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Fig. 1 — Surges when contacts open an inductive load, (a) Test Circuit, (b) Current, 
and voltage vs time after interruption for ideal contact (smooth curve) and typical 
mechanical contact (irregular curves). 
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point of the amount of inductive energy which is available to be dis- 
sipated into the arcing contact circuit, with attendant surges and 
contact erosion. 

1.2 Objectives 

This paper aims to present these relationships in general form, 
with enough experimental data to verify the analysis. Some specific 
fields of relay and circuit design where immediate applications are 
expected, and which were an incentive for this study, will now be 
summarized. 

(z) The general definition of the "ideal" contact surges. It is found 
that all variations are definable in terms of two dimensionless 
ratios: /3 = CR'/(PN 2 , the "load design factor," and G a /Gi or the 
ratio of secondary to primary coil constants (summarized in 
Table I and further described below). The surge equations 
can then be used to predict such variations as peak voltage, 
amplitude and frequency of oscillations, rise time, or degree of 
damping. The analysis also shows that the response is oscillatory 
between two particular values of /3, with overdamped behavior 
above or below these values. 
(ii) An understanding of how the coil's stored inductive energy is 
released. Relationships are found which determine that part 
consumed in internal and external dissipative losses, and that 
part which goes into capacitive energy storage from which it 
may be discharged into erosion-producing arcs. 
(iii) A fuller understanding of that portion of the coil's stored 
inductive energy that can be released during the transient 
period, the major features of which depend on (a) the incremental 
permeance of the static demagnetization curve at the releasing 
point, and (b) the energy replenishment from the battery 
during each arc breakdown interval. This is especially important 
in the study of contact erosion. 
(iv) An explanation of some hitherto "anomalous" results. Among 
these are: (a) A quantitative picture of releasable coil energy 
as influenced by "stop-pin" height; it is shown here to be related 
directly to the degree of coil saturation, in accord with ex- 
perience, rather than to the total stored energy as had previously 
been suggested, (b) A quantitative picture of the oscillations 
which occur with "copper-sleeve" relays. This is found to be 
related to that portion of the main winding which is uncoupled 
to the secondary. 
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Table I — Circuit Constants 

Battery voltage. 

Instantaneous voltage across contacts. 
Steady state current = V„/R w . 

Instantaneous current in primary or secondary circuits respectively. 
Instantaneous charge on capacitor. 
Turns in primary, or secondary, respectively. 
Resistance of primary, or secondary, winding respectively. 
Primary circuit, resistance external to coil (i.e., resistance of RC net- 
work). 

Total resistance of primary circuit = R w + R„. 
Coil constant for primary or secondary circuit 

-= — 5- respectively (mho). 

Eddy current constant of core (for which N« = 1). 

[Note that * = *=+*= = 1+1 ] 

Incremental inductance of primary or secondary 
A r r(P, A'V(P respectively. 

Incremental permeance of coil magnetic circuit at specified magnetiza- 
tion (i.e., corresponding to operated ampere-turns). 
Capacitance across load coil. 
Magnetic flux. 
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(v) The ability to make design estimates for a number of applica- 
tions. Examples are: (a) Number and frequency of breakdowns 
in a "B" type transient, and the energy available for causing 
erosion, (b) Peak value of the final surge, following a "B" type 
transient as for the tail of the irregular curve in Fig. 1. This is 
needed in verifying whether wiring insulation will be damaged 
when coils with large numbers of turns are disconnected. 
(c) Choice of proper test values for use in the production testing 
of relay coils to determine whether their internal insulation 
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breakdown strength is satisfactory. (For this test, which has 
been in use on an empirical basis for many years, the energized 
coil circuit is interrupted by an "ideal" vacuum tube "switch"; 
if the coil's self-induced voltage causes an internal breakdown, 
the coil is rejected as an incipient service failure because of 
faulty insulation or lead configuration.) (d) The ability to 
reinterpret inconclusive data from past contact tests, or to 
better define tests made years ago so that their applicability 
to today's designs in today's circuits may be better understood. 
(vi) Assessment of the mutual impact of contact erosion and relay 
miniaturization. For example, a smaller relay coil usually 
means less energy for the contact to control and therefore a 
smaller contact volume for a given contact life. A smaller 
contact volume means less wear allowance in the relay's ar- 
mature motion, and hence a lower mechanical work requirement, 
with a correspondingly smaller relay magnet. Thus a knowledge 
of the interaction between the relay's inductive energy and 
the resulting contact wear can be an important guide to design 
for minimal relay sizes for a switching system where all the 
relays are of the same basic structure (i.e., the contacts of 
one relay control the winding of another). 
The paper is divided into three major sections: Analysis, Experiment, 
and Discussion. The circuit is first analyzed in general terms, giving 
results which for some cases are considered too complex for further 
treatment in this study. Substitutions are then made which considerably 
simplify the equations and lead to detailed answers for the important 
practical case of the unprotected relay coil. The major experimental 
results and discussion relate to this case. 



II. ANALYSIS 

In this section, we shall (i) develop the circuit equations and apply 
them to the consideration of special situations such as (ii) critical 
damping, (Hi) peak voltage, and (iv) energy relationships. 

2.1 Development of the Circuit Equations 

This section defines a circuit model, presents the loop equations, 
and gives the resulting solution for current in the main circuit. This 
value of current as a function of time is then used to derive the related 
values of charge, contact voltage, and rate of change of voltage. 
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2.1.1 Choice o] Circuit Model 

A simple relay winding shunted by an RC network, being connected 
or disconnected from a battery by means of an ideal contact, was 
chosen for the study. As only the surge on disconnection is studied, 
the influence of wiring can be included in the capacitive term, and so 
the results are generally applicable whether connections between 
contact and relay are very short, as on a printed wiring board, or 
quite long, as through a considerable loop of switchboard cabling. 
Reflections and surge impedance of the wiring are important in con- 
sidering the individual showering arcs, as has been done for example 
in recent work of CJ. W. Mills 10 and J. Pi. Pharney;" but these features 
are outside the scope of the present study. 

After some initial experiments, it was found that an adequate repre- 
sentation of an actual relay control circuit required the recognition 
of: (i) the distributed capacitance of the coil; (u) a coupled secondary 
winding to represent core losses, copper sleeves, etc.; (Hi) an effective 
lumped inductance representing the region of the demagnetization 
curve that is active during the initial surge; and (iv) (under certain 
conditions) a portion of the main winding which is not coupled with 
the core or secondary winding. The resulting study circuit is seen in 
Fig. 2a. For example, (?) it was found that distributed capacitance 
of typical relay coils is in the vicinity of 100-200 pF, which is too large 
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Fig. 2 — Circuits for fully representing a relay as a contact load, (a) Circuit chosen 
for analysis, (b) Transformed circuit. 
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to ignore. Likewise, (ii) it was found that the discharges of typical 
relay structures can be so fast that the core losses (i.e., dissipation 
in an effective secondary winding) must be recognized. The additional 
influence of that portion of the coil winding which is not coupled to 
the "secondary" plays a significant role only when dealing with coils 
having a secondary consisting of a copper sleeve or short-circuited 
turns; it will be examined as a special case of the above circuit in Section 
2.2. Following this general solution, the simpler cases are readily found 
as special cases; e.g., setting the protection network resistance R„ = 
gives the case of an ordinary relay coil with distributed or direct shunt 
capacitance, setting the eddy current coil constant G 2 = gives the 
case of an air core coil, or setting C = °° represents a pure resistive 
shunt. 

The first objective is therefore to evaluate the current, charge, and 
contact voltage in the primary circuit for the period immediately 
after interrupting the circuit shown in Fig. 2a. The following considera- 
tions will be used in the initial analysis: 

(*) The contact is "ideal," i.e., no gaseous or metal vapor break- 
downs occur. This is in fact closely simulated by a mercury 
contact, and even more completely by certain vacuum tubes 
and transistors. 
(ii) Perfect coupling is assumed between primary and secondary. 
(Hi) The initial conditions (just prior to t = 0) are that the steady 
current is /„ = (V„/R,„)* and the capacitor is charged to the 
value — V„ . 
(iv) The flux <f> which is linked to both primary and secondary 
windings is proportional to the algebraic sum of the instan- 
taneous ampere turns of primary and secondary, in accordance 
with the usual magnetic circuit concepts. 
(v) The inductance of the load can be represented by a lumped 
value which effectively simulates the flux-current relation in 
the coil's working region. Use of a lumped-constant term 
simplifies the circuit equations, which otherwise become very 
difficult to solve in simple forms. The meaning of a lumped- 
constant inductance is better understood with the help of 
Fig. 3, which shows the static flux-ampere turn characteristic 
of a typical relay magnet for operating and releasing cycles. 
For the present case, the releasing cycle, or demagnetization 
curve, applies. The instantaneous inductance is the slope of 



For definitions of terms, see Tables I, II, and III. 
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this demagnetization curve in the applicable working range, 
and may be expressed as 

L = <?N a 

where (P = (d$>/dNI) | . v ,„ , N = number of turns in the coil, 
and NI„ = ampere turns at the operated point (i.e., starting 
point of the demagnetization curve on disconnection). (P is 
the differential permeance, or slope of the demagnetization 
curve, over the region where flux is changing during the surge. 
In the case of very fast surges, the flux is substantially un- 
changed during the surge, and the value of (P at NI„ applies, 
as shown at A in Fig. 3. For slower surges, a larger value of (P 
representing an average slope over the working range, as 
indicated by the dotted lines B in the figure, is presumed to 
apply. Values of (P are further discussed in Section 3.1.2. 
The analysis lias been developed around this general circuit. The 
resulting general solution, which involves an appreciable number of 
variables, is relatively complex and its detailed examination is reserved 
for later computer-assisted studies. For the present study, special 
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attention is aimed only at the simpler but important practical case 
of the ordinary unprotected relay, which involves only a simple ca- 
pacitive shunt, a term which is always needed if only to cover the 
inevitable distributed capacitance of the coil itself. In addition, the 
manner in which the general expressions reduce to simpler forms, 
by neglecting appropriate terms, will be noted; conformity to previously 
known solutions for the conventional L, R, C circuit affords at least 
partial confirmation of the general analysis. The special cases chosen 
for detailed consideration are: 

R n = o, G 2 finite (relay coil with capacitive shunt), 
R n finite, G 2 = (air core coil with RC shunt), and 
R n = o, G 2 = (air core coil with capacitive shunt). 
The appropriate expressions for each of these cases are tabulated in 
the accompanying tables. 

The symbols applicable to this model circuit are summarized in 
Table I. The notation using terms G — N 2 /R for coil constant, and 
expressing inductance in terms of turns and permeance (P (L = (PN ), 
has been used because of its convenience and common use in relay 
design. The resulting equations are also considerably simplified by 
their use. 

2.1.2 Equations for Current in Primary 

The equations for the model circuit, Fig. 2a, are derived in Appendix 
A. The procedure was to sum the voltages around primary and secondary 
circuits, respectively, using the flux common to both windings to 
determine the inductive terms, Ni(d$/dt) and N 2 (d$/dt). LaPlace 
Transform techniques were used to represent the initial conditions 
and obtain the solution for current t. in the primary circuit. The three 
possible solutions, oscillatory, critically damped, or overdamped 
[i.e., u 2 a greater than, equal to, or less than b 2 ) are: 



Case I: . ^ Ki e~ bl sin (ut + <p) 

Oscillatory sin <p 



(»l < b 2 ) 



Case II: 

Critically i, = Ki e~ bt [l + (a - b)t] («' - b 2 ) (1) 

Damped 

Case III: . Ki e~'' 1 sigh (fat + v) , ■-> ^ ,2v 

Overdamped ' sinhp 

The values of the constants, of the phase angles, and of related terms, 
are summarized in Tables II and III respectively. The circuit equations 
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also showed that the transformed circuit of Fig. 2b could be used to 
represent all the variations in the primary circuit, where the primed 
values of L and R represent the effective values of primary circuit 
inductance and resistance respectively. This simplification results 
from the term G a <P/C, which is found to represent the dissipative 
losses due to the secondary. (The term L" will be discussed in Section 
2.2.2.) 

2.1.3 Equations for Charge, Current, and Contact Voltage 

The expression for current, equation (1), may be used to derive 
equations for charge and contact voltage as a function of time. The 
charge g, is given by 



and contact voltage by 



h = J h dt (2) 



V e = V + £ + i,R n . (3) 



The resulting general expressions for all these variables are summarized 
in Table IV for the cases: oscillatory, critically damped, and over- 
damped. Equations are also given for two important special cases: 
R„ = 0, and both R n and G 2 = (i.e., capacitive shunt only, or air 
core coil, respectively.) Of course, other special cases can be derived 
by making appropriate substitutions into the general equations. Table V 
similarly summarizes these results for the instant of contact opening, 
i.e., for conditions at t = 0. It is found that, at t = 0, Cases I, II, or 
III all give equal values for each variable, as one would expect. 

The equations in Table IV express the variations in charge, current, 
and contact voltage as a function of time. Of these, it is often of greatest 
interest to consider the contact voltage; a normalized version is plotted 
in Fig. S based on the discussion of peak voltages in Section 2.3. 

2.2 Discussion of Conditions at Critical Dam-ping 

The boundary condition for critical damping, now to be considered, 
helps to explain the way in which the various kinds of surges are in- 
fluenced by the circuit parameters. Comparison of experimentally 
observed conditions at the critical damping point with this analysis 
also provides confirmation for much of the analysis; these measurements 
also provide means for determining effective values for the core loss 
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Table V — Basic Circuit Equations when t = + 



Variable 


General Case 


Case when R n -- 


= 


(a) r/, = 


- CV. 


= -CV„ 


(b) ii = 


KI 


= I, 


(c) V c = 


KI„Rn 


= 


(d) V, = 


A'=/„r CR,* 
C [_ ^i(l +Gi/G a 


,] 


/„ 
c 



term, G 2 , and the effective permeance, (P, of the core, as will be seen 
in Section 3.2. 

2.2.1 General Expressions 

The critically damped condition is ta\ = b 2 , or alternatively: 



CR' 2 _ , 

AL' ~ 1 - 

Substituting values for R' and V from Table I, 

f Ga(P 



(1) 



C \~C~ + Rl ) = 4iV ' (P < 1 + ^t/^i). 
This expression is simplified by making the substitution 

CR 2 



ft = 



<pN'i ' 



yielding a quadratic equation which can be expressed in terms of only 
two variables, /3, and G 2 /G t . The solution is given by either of the 
two equivalent expressions, equations (5a) and (5b), of which the 
former is usually more convenient: 



r; 2 /G\ = $ 03; ± 2) 



or 



ft 



= 9 



v+kA+Wi 



(5a) 



(5b) 



These equations have two solutions, which depend on whether 
G 2 /Gi > ft , or G 2 /Gi < ft . The two solutions are plotted in Fig. 4, 
for the range of these values likely to be encountered in relay switching 
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Fig. 4 — Conditions for critical damping— Go/Gi vs &.* 



practice. Thus, the boundary between oscillatory and overdamped 
conditions is determined by either of two circumstances: 

(7) Losses in secondary circuit are dominant [use + sign in equation 
(5a) or — sign in equation (5b)], for which 



Here, 



f>h, or q f>R 1 - 



P 2 = ti(p\ + 2) 



(6a) 



or 



0. = 2 



■+*-'+£/ 



(6b) 



and there are positive values of f/ 2 /(7, for all positive values of ft . 



* The subscripts "1" or "w", applied to G and /3 in various figures, are used to call 
attention to the applicability of the information; e.g., ft indicates that the general 
analysis (including RC network) applies; ft,, indicates the analysis specifically applied 
to the unprotected case, where R n = 0. 
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(«) Losses in primary circuit arc dominant [use — sign in equation 
(5a) or -f sign in equation (ob)], for which 



Here, 



7T,< 



( ll _ , 



, or R x > 



r/ 2 (P 
C ' 



g = WW - 2) 



or 



ft 



1 + ^ + ( 1 + § 



(7a) 



(7b) 



and meaningful (i.e., positive) values of G 2 /G x only exist when 

ft > 4. 
Note that, when G 2 = 0, i.e., there are no eddy current losses, then 
from equation (5b), 

ft = 4 

as is well known for the conventional L, R, C circuit. 
The condition 



G, 



= ft l or 



C ~ K] 



is seen to be incompatible with critical damping. It is indicated by the 
45-degree line in Fig. 4 which in fact corresponds to the point of mini- 
mum damping. The character of the contact voltage surges in the 
various regions is sketched on the figure, summarized as follows. For 
a given value of G 2 /G\ , the response will be oscillatory for values of 
fix lying between the two curves, with minimum damping when ft = 
G 2 /G, . When 0, equals or exceeds the value delineated by the right- 
hand curve, the response will be critically damped or overdamped 
with a gradual asymptotic rise to Y . When /3, equals or is less than 
the value delineated by the left-hand curve, the response will be critically 
damped or overdamped with a single overshoot followed by an asymp- 
totic decay to V„ . 

Since the principal variation in /3, is usually the result of variations 
in C, it is also useful to replot Fig. 4 in the form G 2 /G, vs C, for a range 
of values of {N x /R i )'(9. Such a family of curves is given in Fig. ">. 
It is sometimes helpful to analyze experimental data by plotting the 
points on such a graph. 
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Fig. 5 — Conditions for critical damping — G-/Gi va C. 

The data points seen in Figs. 4 and 5 will be explained in the later 
discussion of experimental results, Section 3.2. 

The critically damped circuit equations, Case II in Table IV, may 
now be extended to include expressions for the two-valued equalities 
of equation (5), which were not previously obvious. For this purpose 
it is convenient to make available the related values of some of the 
terms appearing in these equations. In the present study we consider 
only the simpler cases for R n = 0, or both R n and G 2 = 0, which are 
summarized in Table VI. Using these equalities, the critically damped 
circuit equations of Table IV can be simplified for easier handling, 
as given for reference purposes in Table VII. 

The predicted boundaries between overdamped and oscillatory 
conditions shown in Figs. 4 and 5 afford a convenient means for experi- 
mentally verifying the foregoing circuit analysis. A test relay is selected 
which has parallel-wound primary and secondary coils. The primary 
is bridged by a variable capacitor, C, while the secondary is bridged 
by a variable resistor, R!, . For a chosen setting of C, the contact voltage 
surge V e is observed by oscilloscope, and R' 2 is varied until the surge 
appears to be critically damped. This process may be repeated for 
any desired number of settings for C, which must include corrections 
to include any distributed capacitance in the winding and wiring, 
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as is discussed in Section 3.1. These data may then be analyzed to give 
estimates for G, and <P for the particular structure being studied; 
and also to compare the measured variation of 6*2/6', vs /3 against 
the predicted values, as is discussed in Sections 3.1.3 and 3.2. 

In accord with relay design practice, assume that the effective 
secondary coil constant, 6' 2 , is the sum of two separate terms, 

g. = v 



H, + R' 2 



and G f , representing the actual secondary winding, and the single- 
turn eddy current coupling term for the core, respectively. An ap- 
proximate estimate of G, may then be obtained from an examination 
of those data points for which G, has quite small values. This procedure 
is helped by the fact that, at small values of G 2 /Gj and /3, the critical 
damping equations (6) may be further simplified. By series expansion 
of equation (Ob), neglecting terms of higher than second order, it 
reduces to 

or, from Table I, 

G 2 = G, + G, = 22V, VCAP. (9) 

Thus, a plot of G, vs C*, for small values of 6' 2 , should yield a straight 
line with a negative intercept equal to G, and a slope of 2iV 1 <P~*, from 
which (P may also be determined. This procedure is indicated in the 
experiment summarized in Table VIII, and further discussed in Section 
3.1, with determination of G, shown in Fig. 17. In this example, G e 
was found to equal about 220 mho. Using this figure to obtain an ap- 
proximate value for 6' 2 , the effective value of permeance (P (at least 
for small values of 6' 2 ) may then be determined using equation (9), 
as is also shown in Fig. 17. 

2.2.2 Peak Voltages and Oscillations due to Imperfect Coupling Between 
Primary and Secondary 

The foregoing analysis is found to afford an accurate picture of the 
transients upon disconnection, so long as the windings are well coupled, 
as is shown in Section III. Such coupling is best obtained by the use 
of parallel windings, wherein each turn of one coil is coupled to the 
same flux as the corresponding turn of the other coil. However, when 
the coils are wound concentrically, there is a departure from this ideal 
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situation due to independent transient currents which arise within 
the uncoupled portion of the main winding. An additional oscillatory 
condition develops in such cases, which is particularly noticeable for 
relays with copper sleeves, whose core losses tend to be very high; 
this condition becomes more and more noticeable when the circuit 
approaches critical damping or is well into the overdamped core-loss 
region. This behavior may be explained by recognizing the role of the 
uncoupled part of the primary winding, now to be discussed. 

Recall that the core loss is governed by a term G 2 (P/C representing 
the resistive losses in the secondary, and note that this term grows 
very large as C becomes small. Thus, in the critical damping region, 
this large term tends to block all flow of current in the primary in- 
ductance. However, that portion of the primary inductance which is 
not coupled to the secondary, i.e., which does not share the same common 
flux, will not experience opposition to current flow from the term 
G 2 (P/C (because, for it, (7 2 = 0). The uncoupled part of the winding 
provides a parallel path tending to bypass the main winding in series 
with its core-loss term, as shown by the dotted lines in Fig. 2b. Its 
inductance may be much less than that of the coupled part of the 
winding, but still capable of generating appreciable superimposed 
surges into the rest of the circuit. The equivalent circuit for this con- 
dition reverts to the simple case for which only coil losses apply (i.e., 
no secondary winding) but for which the inductance is represented 
by a new term 

L" = <P"Nl (10) 

where (P" represents the effective differential permeance of that part of 
the coil which does not couple to the secondary. It may be found 
experimentally, or be approximated by rule of thumb, as discussed 
in Section III. The resulting oscillations, as to frequency, peak voltage, 
etc., may be predicted through the use of the equations in Table II, 
for the case that G 2 = 0, and using L" to represent the inductance. 
Peak voltage effects for this case may be read from Fig. 8, where fi 
is estimated by using (?" for permeance. The oscillation frequency 
is approximated from the value of o> given in Table II: 

_ w « 1 

/o = tc = 2*N l V¥ 7 C' (11) 

Conversely, when the oscillation frequency is measured, corresponding 
to a given set of circuit conditions (i.e., 0, G, etc.), the value of (P" 
can be evaluated. Measurements on various relays are shown in Section 
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III to correlate well with this relation, where (P" is found to have a 
value of about 0.6 X 10" 8 Wb/NI for various structural types typical 
of telephone switching. 

2.3 Determination of Peak Values of Contact Voltage, and Related Values 
of Charge and Current 

The peak values of the time-dependent variables q t , ii , and V e 
may be determined by differentiating the corresponding equations 
as given in Table IV with respect to t, setting the result equal to 0, 
solving for t under these conditions (designated t'), and then evaluating 
each of the original equations when t is replaced by t' . These steps 
together with the resulting expressions for t' and for </i , i t , or V e are 
briefly outlined in Appendix B. The derivations apply to the previously 
noted general case which includes an RC network shunted across the 
load, but the resulting general expressions are too involved for a full 
treatment in the present report. When the network resistance R n is 
set equal to zero, however, representing the simpler unprotected case, 
the resulting expressions are relatively convenient to work with. These 
expressions, which are also developed in Appendix B, are summarized 
in Table IX. The peak voltage relations hold greatest interest, and 
will now be considered in more detail. 

As seen in Appendix B, and Table IX, the peak voltage expressions 
for an unprotected load (R n = 0) are: 

Contact 
Condition Time to Peak, t' Peak Voltage, V p , 

ir.z=£ Fp . F .[ I + <c .„ p (_= i ^)]. cud 

II = (i,Ml(l + A)' = V.[l + p-Je-"" 1 -' ]. (13) 



i" --S =4 1 + ft M-u^)J- (14) 

These equations express the peak voltage in terms of the variables 
/3„ , 0, and tp. Actually, both 6 and <p may also be expressed in terms 
of p m and G 2 /C w (see Table III). The equations are given in the form 
shown for the sake of compactness. However, when computations 
are made, the more detailed values for these terms can be substituted 
from Table III. (The subscript w is used to indicate that the value 
of R in this case is only the winding resistance, i.e., R„ = 0). 

For this case of unprotected load, equations (12), (13), and (14) 
may be used to plot the time-to-peak, and the peak voltage. It is 
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Fig. 6 — Time, I', to reach peak voltage vs load design factor, 0. 

convenient to plot these in ratio form, using the ratio of peak time to 
capacitive time constant, t'/R w C, and the ratio of peak voltage to 
battery voltage, VJV B , respectively. These results are given in Figs. 6 
and 7, as a function of fi. w , for various values of G 2 /G n . They are 
found to be in good agreement with actual experience, as will be shown 
in Section 3.3, where measurements and theory are compared for both 
t' and V p . It is seen that the ratio G 2 /G a has a pronounced influence 




Fig. 7 — Peak voltage, V p , as a function of load design factor, 0. 
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on the peak voltage across the opening contacts, and that even very 
small values of this ratio can lower the peak voltage and pulse time by 
several fold compared to the air core (i.e., nondissipative) case. This 
is particularly noticeable when small values of W (usually corresponding 
to small values of shunt capacitance) are involved. 

A summary of the equations for the peak-related values of g, , i x , 
and V c is also given in Table IX. 

The relations determining these peak values may be used to restate 
the circuit equations in ratio form, if desired. As the contact voltage 
expressed in this form is of particular interest in later discussions, 
its value (when R n = 0) is given below as 



V r - Vo 
V — V 



pxp 



(-Sfo-4 



, m f (1 - r))sin {(tt - ^)r + ^ + 6 



sin 6 



(15) 



Equation (15) expresses the ratio of instantaneous to peak coil voltage 
as a function of the ratio of instantaneous to peak time, t, following 
disconnection. A plot of this expression for two values of W is given 
in Fig. 8. 

2.4. Energy Losses During the Surge 

The surges just discussed can lead to breakdowns in wiring, in 
apparatus, or between the opening contacts. The breakdown itself 
is usually a very rapid event resulting from the discharge of the circuit 
capacitance through the low-impedance breakdown path of electrodes 
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Fig. 8 — Normalized plot of coil surge voltage (V c — V„) vs time. 
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and battery. The energy transferred from inductance to capacitance 
is thus the preponderant source of breakdown damage, and expressions 
are therefore needed to describe the manner in which the inductive 
energy is dissipated or transferred to the capacitor. In this section, a 
brief discussion is given of the breakdown process, followed by the 
derivation of energy expressions based on the foregoing circuit equa- 
tions. 

2.4.1 Breakdowns 

The earlier expressions for contact voltage are of particular interest 
when breakdowns occur, because it is the value of voltage — matched 
against the dielectric strength of the wiring, of the connected apparatus, 
or of the contact itself — which determines when and how often break- 
downs will occur. These breakdown possibilities are indicated in the 
simplified view of Fig. 9, which shows a typical voltage surge vs time 
after disconnection, and also two forms of dielectric characteristics 
representing the breakdown strength of portions of the circuit. The 
horizontal line, designated A, represents, for example, the breakdown 
characteristic of a pair of wires or apparatus terminals with constant 
spacing. The line of rising slope, designated B, represents the break- 




Fig. 9 — Interaction between inductive surge and contact's dielectric characteris- 
tic (simplified view). 
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down characteristic of separating contacts. Its characteristic is well 
approximated by the equation 

V B = Y 00 + kut (16) 

(see for example Ref. 9) 

where V 00 = minimum gas breakdown (usually the Paschen breakdown 
point which for air is about 320 volts), 
k = slope constant (for air it is approximately 7 X 10 6 V/m), 

and 
u = contact separation velocity. 
There will be a gas breakdown whenever the coil's voltage characteristic, 
plus the battery voltage, exceeds the electrode separation characteristic, 
as is shown in the figure for the case of separating contacts. Each 
breakdown is commonly followed by an arc resulting from the discharge 
of the local capacitance, which was charged to a voltage (V B — V ), 
through the low-impedance local wiring, usually causing a partial 
capacitive discharge and an extinction of the arc due to reverse oscil- 
lations. Following this, a new charging cycle begins which is similar 
to the earlier surge except rising toward a somewhat lower peak voltage 
representing the energy remaining in the load after the first breakdown. 
This procedure will continue repetitively until the residual load energy 
is insufficient to develop a surge whose peak value can reach the elec- 
trodes' dielectric breakdown. 

For the constantly-spaced electrodes, each disconnection will lead 
to a group of multiple breakdowns which induce radiation in the 
wiring and a slow erosion of their surfaces. For example, fine-wire 
coils with too-close leadouts can eventually develop open-circuits as 
a result of this slow attrition. 

Similarly, separating contacts experience an erosion of their surfaces. 
In the case of relay contacts, it is desirable to correlate the amount 
of erosion with the circuit parameters in order to determine either 
(?) the required contact volume to permit a given number of operations 
before failure, or (u) the number of operations for a given contact 
volume. The remainder of this section will indicate one approach 
toward the correlation of contact erosion with the conditions of circuit 
and contact. 

Many contact erosion studies have shown that the amount of metal 
lost is proportional to the energy dissipated between the electrodes 
during the discharge. 7,9 For the present case, a method of estimating 
this energy for the multiple breakdowns of the "B" type transient is 
therefore needed, and two alternate procedures are conceptually 
possible. The two approaches may be more clearly understood by a 
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reexamination of the oscilloscopic traces of coil current and contact 
voltage immediately following the contact opening, in Fig. 1. Each of 
the multiple surges for the mechanical contact is seen to be the early 
stage of an ideal contact surge which is interrupted when it reaches 
the breakdown characteristic of the contact gap, as is also indicated 
in Fig. 9. The breakdown is due to the voltage to which the capacitor 
is charged, plus the battery voltage, and when the breakdown does 
occur, it represents the discharge of this capacitance through the 
circuit comprising capacitance, contact, battery, and the low impedance 
of their wiring. This is a very rapid discharge, occurring in times that 
are negligible compared to the time required for the original surge to 
reach its breakdown value. The discharge is extinguished either by 
exhausting the capacitive stored energy \C(V B — V„) 2 , or by reversals 
in voltage and current during the very high frequency breakdown 
discharge. Following extinction, a new surge develops as before. The 
energy dissipated in these multiple arcs may be approached in two 
ways: (i) by summing the individual energy terms (each of the form 
&E = j l V A I A dt), where V A is the characteristic voltage across the 
arc for the particular contact metal (typically about 14 volts) and I A 
is the current through the arc, or (ii) by the total of the energy in- 
crements which are transferred from the coil to the capacitor to the 
arc (which should be related to the total initial inductive energy). 
Though the determination of V A , I A and t are subject to many dif- 
ficulties of estimation, efforts to treat this problem with computer 
techniques have shown much promise. 10 The alternative procedure, 
namely to find the inductive energy which can be released into the 
arc, via the capacitor, will be considered here. Contact erosion should 
be strongly related to this energy, which will be discussed for the ideal 
inductive case in Section 2.4.2, and for the usual showering case in 
Section 2.4.3. 

A more rigorous treatment must include the additional energy 
supplied from the battery during each breakdown arc, and an intro- 
ductory discussion is given in Section 2.4.4. 

2.4.2 Energy Dissipation Following Disconnection — Ideal Contact 

First, let us find the energy consumed during the inductive surge, 
as a function of time. For simplicity, the analysis given here will con- 
sider the case that R n = 0, i.e., the unprotected load. Based on the 
equivalent circuit of Fig. 2b, and equation (lb) Col. 2 of Table IV, 
the energy dissipated in the circuit resistance is seen to be 
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E reB = f i 2 R' dt = 4^- f e" 2 "' sin 2 U + <p) dt. 
J„ sin (p J 



(17) 



Upon integrating between the limits indicated, this internally dis- 
sipated energy is found to be 

T'T 2 

E res = 7T^~ [U - C OS 6 COS (2<f + e)\ 

z sin ^ 

- e -z "{l - cos e cos [2(«* + <p) + 5]}]- (18) 

A closer inspection of equation (18), and further rearrangement of 
terms based on the relations of Tables II and IV, gives the alternative 
expression 

#re. -^(V+ajG. + M 

_ £ - 26 ,r sin 2 (ut + <p) + sin 2 (ut + g + 0) 1 .-„» 
L sin 2 J 

= A'. - #„ (20) 

where ^„ = initially available energy in coil and capacitance 

= wn + hcvi , (2D 

E u = undissipated energy, i.e., that portion of initially stored 
energy, which, at time t, has not yet been dissipated 

= E ind , energy still stored in inductance, L', 

+ E cav , energy so far transferred to capacitor, C. 

These expressions thus offer the means for quantitatively estimating 
the way in which the initially stored energy is either internally dis- 
sipated (in the resistive term R') or made available for release outside 
the system (i.e., stored in the capacitor). 

The extent to which energy may be released outside the coil system, 
as for example into a contact arc, is perhaps more easily pictured by 
further reference to Fig. 2, recalling that, for the present discussion, 
R n = 0. At the time when the capacitor C is charged to a value capable 
of producing a dielectric breakdown at a contact voltage, V B , it possesses 
an energy content hC(V B — V n ) 2 which can be rapidly dissipated in 
the low-impedance contact circuit during a breakdown. The maximum 
capacitive energy which can be so discharged occurs at the time t' 
and at the voltage V p . t' is the time when the main circuit current 
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equals zero; at this instant all the undissipated inductive energy has 
been transferred into the capacitor. The energy in the system at this 
point, E' u = E caf> , is therefore a useful measure of the maximum portion 
of the coil energy which can be released to do external damage. By 
substituting the value of t' from equation (12) into the second right- 
hand term of equation (19), this is seen to give 

2(t - <p)' 



E 



cftp 

E 



EL 

E n 



exp 



tan 9 



(22) 



1 + G 2 /G w + ft. 

Equation (22) states the maximum fraction of the initial releasable 
inductive energy that would be capable of doing external damage. 
Since <p and are themselves (from Table III) expressible solely in 
terms of the circuit and relay design parameters j8 w and G 2 /G w , this 
energy factor is similarly dependent on only these two variables. The 
result for a wide practical range of relay and circuit variables is given 
by the family of curves plotted in Fig. 10. It is interesting that these 
curves indicate a central range of values for /3„, where the energy that 
can be externally released is a maximum, with lesser energy values 
for larger or smaller values of /3 1( , . This seems physically reasonable 
because small values of /3„ correspond to small values of C and rapid 
pulses, while large values of /3„, correspond to slow pulses. Thus, for 
the short pulses, more energy is forced into eddy-current dissipation, 
while for the slow pulses there is ample time to dissipate energy in 




0w= 



C R x 
<PN. 



Fig. 10 — Fraction of releasable energy, E eap /E , available as potential source of 
arcing damage vs /3. 
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the coil resistance before the capacitance is charged; for intermediate 
values of (3„ , both dissipation terms are relatively small. 

To summarize this concept of maximum releasable energy, consider 
an illustrative case for which the following circuit and relay design 
conditions are given: 

Circuit: V = 50 volts, C = 10" 9 F, I„ = 0.100 ampere. 
Relay Coil: JV, = 3100 turns, R u = 500 ohms, 

(P = 2.5 X 10" 8 Wb/NI, G 2 = 100 mho. 
Determination of parameters and associated energy: 

_ 1Q-'" X 250,000 _ : , 
Pw 2.5 X 10 -8 X 10 a " 

G w = i^- = 20,000 mho; G 2 /G„ = 0.005. 

E n = — * IU (3100 X 0.1)- 

= 1.25 X 10" 3 J, initially available. 

From the chart: E cap /E = 0.74, or maximum releasable energy = 0.92 
X 10 -3 joule. Since contact erosion is believed proportional to the energy 
delivered to the arc, this number is a useful starting point in establishing 
the system's erosion-inducing energy. For arcing contacts, it may be 
materially improved, it is believed, by an upward correction to the 
initial energy E„ as will be discussed in Section 2.4.4. 

2.4.3 Energy Dissipation During "B" Type Transient 

The preceding results are directly applicable to the repetitive break- 
downs that occur between ordinary relay contacts that are separating 
in an air environment. The breakdowns occur when the surge voltage 
transient exceeds the dielectric breakdown characteristic of the contact 
gap, as previously discussed in Section 2.4.1 and illustrated in Fig. 9. 
Each breakdown occurs during a surge of the form of equation (15), 
which rises to intersect the contact characteristic given by equation 
(16) as illustrated in the figure. 

2.4.3.1 Steps During (he Breakdown Process. The energy losses during 
this series of breakdowns will thus be somewhat different from those 
in the idealized case discussed in Section 2.4.2 because they occur at 
an earlier stage of the surge. This effect may be visualized by con- 
sidering the following simplified statement of representative sequential 
steps in the total breakdown process: 
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(?') The surge rises to a value of V B given by equation (16). An 
amount of energy §C(V B — V„) 2 is released into the arc circuit 
while concurrently an amount of energy given by equation 
(19) is dissipated internally. The energy still remaining in the 
inductance to serve as a source for the next surge is E ind and it 
comprises the initial energy of the previous surge minus the 
capacitive and resistive terms just mentioned. (A later and 
more rigorous study must also include the capacitive energy 
that is not dissipated during the surge, and so temporarily 
returned to the system, and also the energy replenishment 
from the battery.) 
(ii) A new surge now develops which will rise toward a lower peak 
value defined by the inductive energy content after the previous 
surge, and a new breakdown occurs when a higher value of V B 
is reached, as determined by the intersection of the coil surge 
with the contact's dielectric characteristic, equation (16). 
From the earlier work, 9 or from an inspection of Fig. 9, each 
new breakdown voltage is approximated by 

^-M» + T5&J-Mr=wJ » 

These approximate expressions result from assuming (a) that 
the relevant part of the surge is closely represented by using 
its initial slope, i.e., V/t = IJC when t = (from Table V, 
Case 1(d)), and (b) that the average initial current from the 
first to the last of the series of surges is given by I„/2. 
(Hi) Similar steps continue until the amount of energy that must 
be imparted to the capacitor in order to produce a breakdown 
exceeds the energy remaining in the inductor. Further break- 
downs cannot occur, and any remaining energy is dissipated 
in the coil's internal resistances. 

2.4.3.2 Energy Relations During the Surge. An approximate method 
of estimating the energy that will be available as a driving source for 
arc damage will now be given. (One method would be to use a series 
of linear approximations for energy similar to that for V/t in item ii 
of Section 2.4.3.1. However, the procedure now to be given offers some 
simplicity, more rigor, and more insight into the mechanism of energy 
loss.) The method chosen requires a more detailed picture of the energy 
dissipation during the first stages of the surge, and a statement of the 
necessary energy content at each successively increasing breakdown 
voltage, now to be examined. 
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Equation (19) makes it possible to compare the energy lost in re- 
sistive elements, the energy still remaining in the inductor, and that 
stored in the capacitor, as a function of time. Thus, from equation (19), 

E 

-J*' (24) 



E^ _ _ Em 
E„ E n 



It is readily verified that the first and second (sin) 2 terms in equation 
(19) define the inductive and capacitive energies, respectively. Nor- 
malizing these expressions, as was done for V e in Section 2.3, we then 
write 



E 

E^ 
E„ 



2(tt - <p) 

exp \ — \ f 1 t 

tan 



sin 2 {(» - <p)t + tp\ 



(1 + G a /G w + ftO sin 
exp (- 2( [ a ~/ ) r) sin 2 {( T - v )r + <p + d\ 



(25) 



(26) 



(1 + G 2 /G w + (8.) sin" 5 

Equations (24), (25), and (26) are plotted in Fig. 11 for the case that 
G 2 /G w = 0.005, and two different values of 0„ (= 10" 4 and 10" 3 ). 
These curves show that the total losses during a discharge series may 
be divided between resistive and capacitive losses in quite different 
ratios, depending on the time at which the breakdown occurs, and 
also on the values of /3 U , and G a /G a . Thus, arc damage at the contacts, 
which comes only from the capacitive energy term, E QaD , will be a 
larger or a lesser fraction of the total initial inductive energy, depending 
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Fig. 11 — Dynamic energy relationships during inductive surge. (G>/G w = 0.00. r )). 
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on its ratio to the energy which is concurrently dissipated resistively, 
E Tea . By means of the relations given here, these various forms of 
energy loss can now be calculated, based on the known properties of 
the control contact, the relay load, and the circuit between them. 
In illustration, it is evident from Fig. 11 that the resistive losses com- 
prise a much higher fraction of the total initial energy for the lower 
value of /3,„(10 -4 ). Accordingly, contact wear should be less for this 
case than for the case when /3„ = 10" . 

2.4.3.3 Illustration. The following example illustrates the steps for 
making a simple estimate of energy into the arc for a particular case. 
The needed starting data are as follows: 

|~Control contact: separation velocity, u; V B vs gap as given 
Given: by equation (16). 

.Circuit and load: R„ , N r , C, G 2 , V , NI , (P, (R n = 0). 
Estimation of Energy: 

(i) From above data, find G„ = (N 2 JR W ), G 2 /G w , and /3„. . 
(w) Based on (3 W and G 2 /G w , find V p /V from Fig. 7, which in turn 

determines V p . 
(iii) The values of u, k, I, and C give a ratio which determines each 
successive new breakdown voltage according to equation (23). 
This in turn determines the capacitive energy which must be 
drawn from the inductive source, for each successive breakdown. 
(iv) Reference to the normalized surge curves (Fig. 8) gives the 
normalized time for successive surge breakdowns, from which 
the values of E cat> , E TCS , E ind , and residual E ind can be found 
(Fig. 11). 
(v) This procedure is continued iteratively until the residual 
inductive energy is less than the required capacitive energy 
for the next sequential breakdown. At this point, no further 
breakdowns can occur. The cumulative capacitive energy ratio 
for all the breakdowns to this point then determines that 
portion of the initial energy which can produce contact arc 
damage. 
These steps are of course conveniently performed by computer, but 
are not too complex for approximation using a slide rule. In one trial 
example: a relay with values of G 2 /G w = 0.005 and /3„, = 10 -3 , and 
other constants as in the example of Section 2.4.2, a slide-rule cal- 
culation indicated a series of 10 breakdowns, which checked with 
experiment, and gave at the time of the last breakdown: 

Cumulative %* = 0.34, 
E„ 
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Residual ^ = 0.11, 
hj 

E 

Cumulative —~ = 0.55. 
E a 

Thus, the indicated energy which is capable of producing contact 
damage in this case was estimated to be 

E eav = 0.34 X E 

= 0.34 X 1.25 X 10" 3 

= 0.425 X 10" 3 joule. 

In the specific example in Section 2.4.2, it had been estimated that the 
maximum releasable energy that could be expended in arcing was 
0.74 E . However, in this extension of that case, where multiple break- 
downs occur, with attendant changes in the ratio of resistive to ca- 
pacitive losses, and with some energy still remaining in the system 
after the last breakdown, the energy releasable into the arc circuit 
is found to be. only 0.34 E„ , i.e. about 46 percent of the hypothetical 
maximum. In general, it is expected that the energy releasable to 
cause arcing erosion will be considerably less than the maximum energy 
as given by Fig. 10, though roughly proportional to it. 

As already discussed in Section 2.4.2, the above procedure ignores 
a replenishment of energy into the system, during each brief arcing 
interval, due to the closed-circuit paths for battery current into the 
capacitor and into the inductor. Correction for these effects, discussed 
in the next section, will modify the results just obtained; it is expected 
that an upward revision in E a , which recognizes the circuit current 
as well as the initial inductive energy, will result. The stepwise analysis 
just given is therefore an approximate treatment. Its special value is 
the insight it offers into the apportionment of energy available at the 
onset of each surge between internal dissipative losses, and between 
capacitive and inductive portions, which are the potential sources 
for arcing and associated contact erosion. 

2.4.4 Modifications in Energy Term, E„ , Due to Replenishment of 
Inductive Energy During Showering Arcs 

When the capacitive energy is released into the arc, at breakdown, 
oscillations in the arc circuit tend to extinguish it before the full energy 
is released — either because contact voltage drops below the arc voltage 
or current falls below the minimum arc sustaining current. The un- 
released energy remains in the system; some will be dissipated in the 
circuit resistance and some will be transferred to inductance and back 
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to capacitor, and so will contribute to later breakdowns. The total 
of such available energy has an upper limit already given by E' u in 
equation (22) , of the form 

E' u = E j(t) 

where E = (P(NI ) 2 /2 and (P represents the differential permeance 
at NI„ . But when additional energy is pumped back into the system, 
from the battery, during each arc, the term E„ must be replaced by 
(E + Ei) where E t is this incremental battery energy. Information 
about, the term E { is therefore essential to a full understanding of 
energy entering into the arcing and eroding process. 

A first indication of the approximate values for E { comes from 
empirical erosion studies such as those discussed in Ref. 9. In this 
work it was found that E t = kl where I is the steady state circuit 
current and k is a term related to (i) arc voltage (ii) total time of all 
the arcs, (Hi) magnetic constants of the circuit, etc. For a particular 
group of telephone relays it has the approximate value k = 0.1 (in 
volt-seconds) . 

Further treatment of this correction term, which can be a major 
factor in some cases, is beyond the scope of the present paper. But 
because of its probable significance in many practical cases, its existence 
should be recognized. An extensive study of the importance of this 
battery-contributed energy replenishment, now nearly complete, by 
Pharney 11 , confirms the importance of a corrective current term as 
indicated above. 

III. EXPERIMENT 

The foregoing analysis has been found to agree well with experiment 
for a variety of applications problems. Some of these experiments 
will now be summarized, indicating (i) experimental determinations 
for some of the important circuit constants, (ii) confirmation of various 
critical damping relations, and (Hi) verification of the peak time and 
voltage relations. While the data apply specifically to the case of an 
unprotected load (R„ = 0), the methods seem equally applicable to 
future studies of the more complicated general case. 

3.1 Evaluation of Circuit Parameters Which Describe the Relay Structure 

The principal variables to be considered are the effective distributed 
capacitance of the coil, its effective inductance, and the core-loss term. 
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3.1.1 Distributed Capacitance, C d 

When one measures the voltage surges which develop upon dis- 
connecting a coil in a circuit with very short wiring and no other ap- 
parent capacitance, it is found that the circuit behavior is as though 
a capacitance of the order of 50 to 200 pF were bridged across the 
coil. This is called the coil's effective distributed capacitance, C d ; 
usually its value is fairly constant, for a particular basic structure 
and form of winding, regardless of number of turns, size wire, etc. 
This value should be known with fair accuracy in order to predict the 
relay's surge behavior in an actual circuit. 

The initial slope of the surge voltage curve may be used to find this 
value of apparent capacitance for any particular coil structure. Note 
in Table V that, for the case that R n = 0, the value of V c at time t = 0, is 



dV r 
dt 



Thus, to estimate the value of C, one may open the coil circuit, prefer- 
ably by means of an ideal contact, and note the initial slope of the 
voltage surge, which can be quite easily measured with an oscilloscope. 
Then the indicated value of shunt capacitance is 



C = - 1 .- 



(27) 



This technique can be used to measure any given structure, when all 
external capacitance is carefully minimized. To avoid possible errors 
from use of a single measurement, it is sometimes worthwhile to repeat 
the above measurement after adding several different values of external 
shunt capacitance. Figure 12 demonstrates this procedure in a typical 
measurement made on a particular form of relay structure. The seven 
approximately straight lines on the oscillogram (Fig. 12a) are the 
initial slopes when the bridged capacitance was successively changed 
from to 50, 100, 200, 400, 750, and 1000 pF. The upper trace is a 
measurement of current taken simultaneously, to give the proper 
value of /„ for substitution in equation (27). After making appropriate 
corrections for capacitance of the measuring leads, the apparent ca- 
pacitance was calculated and plotted against the actual values used 
in the experiment. This plot is seen in Fig. 12b. The intercept on the 
apparent -capacitance scale corresponding to C extern:il = is seen to 
be 140 pF, which is taken as the effective value of C d for this particular 
test coil. 
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Fig. 12 — Experimental determination of distributed capacitance, CV (a) Measure- 
ment Of V c and I a . (b) Plot Of Cappurcnt VS Cextennl- 



Values for C d which were measured for a number of different structural 
forms of relays are summarized below, together with brief data on 
relay construction. 

Measured Values of Cd for Certain Relay Coils 









Coil Description 








Code 






















Designa- 






Relative 










tion 


Ni 


Ru, 


Volume 


Shape 


Form of Structure 


c d 


286 


2,740 


288 


large 


1/2 full 


wire spring 




120 


280 


26,550 


12,450 


large 


full 


polarized 




125 


AF33 


11,850 


950 


medium 


full 


wire spring 




140 


AJ5 


19,400 


2,500 


medium 


full 


wire spring 




125 


AJ7 


34,900 


9 , 100 


medium 


full 


wire spring 




110 


AJ12 


5,150 


700 


medium 


short 


wire spring 




110 


AJ57 


P 3,235 


200 


medium 


concentric 


wire spring 




140 




S 3,235 


200 


— 


— 


— 




— 


AJ503 


2,110 


270 


medium 


hollow 


wire spring 




200 


AG55 


P 6,525 


1,000 


medium 


concentric 


wire spring 




180 




S 2,260 


42 


— 


— 


— 




— 


MA24 


P 2,400 


200 


small 


full 


flat spring 




115 




S 2,400 


200 


— 


— 


— 




— 


MB9 


P 4,200 


440 


small 


full 


flat spring 




170 




S 3,500 


350 


— 


— 


— 




— 


295A 


5,180 


525 


medium 


full 


6 large mag. 
contacts 


reed 


195 


311A 


3,950 


550 


small 


full 


3 small mag. 
contacts 


reed 


210 


327B 


7,740 


950 


small 


full 


1 small mag. 
contact 


reed 


140 


327D 


16,550 


2,000 


small 


full 


1 small mag. 
contact 


reed 


150 


Miniature 4,300 


1,560 


very small 


full 


armature 




80 


Structure 
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3.1.2 Effective Inductance 

The present study has postulated a circuit model which comprises 
lumped parameters, including the inductance. The inductance of an 
air core coil should rigorously satisfy this assumption. However, for 
coils using ferromagnetic materials, the flux linkages are a nonlinear 
function of the coil current. It is therefore necessary to consider the 
magnetization curve of the load circuit and to determine the effective 
working region for the particular load structure and circuit condition. 
This may be clarified by further considering the magnetization char- 
acteristic for a typical load structure shown in Fig. 3. The dotted 
curve shows a representative path that may have been followed by 
the prior operation of the unit, and point A indicates the starting 
condition of the inductive load at the time it is disconnected. The 
heavy curve is the demagnetization characteristic of the structure, 
which is the region of interest for the present discussion. Now, if the 
load is connected in a fast circuit (as for example with only a small 
capacitance) the circuit equations determine that the transient re- 
sponses will be so fast that the core flux cannot change appreciably 
(being retarded by the core's eddy currents) in the time available 
during the early stages of the transient. In this case, the initial slope 
of the demagnetization curve (i.e., slope at Point A) should define the 
inductance quite accurately, and one may represent it as L = (?N 2 
where (P represents the differential permeance or slope of the demag- 
netization curve at the operated point, i.e., 

d$ 
(P. = 



dNI 



However, if the load is in a slow circuit (for example with a large 
capacitance), there can be time for an appreciable amount of demag- 
netization. The inductance then grows larger as the transient progresses 
(i.e., (P increases as the $-NI characteristic is traversed downward). 
The assumed lumped inductance no longer rigorously applies. 

In the present study, for the sake of analytical simplicity, this 
variation will be either ignored, or be treated by assuming a progres- 
sively larger value of (P as conditions combine to give slower transients. 
Based on the lumped-constant approach, (P would be estimated by 
an averaging procedure over the portion of the demagnetization curve 
which is traversed during the transient. While it is believed that (P RT 
will largely depend on the ratio of the peak transient time t'/RC to 
the core-loss transient time-constant, (? 2 (P, the subject is not treated 
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further here except as an averaged value of (P is determined empirically. 
An averaged value of (P is indicated in concept by the slope of dotted 
lines such as B in Fig. 3. Use of (P„ usually gives good agreement with 
the analysis over several decades of variation in /3. 

(P may be determined by magnetization data on the particular 
structure, or by various interpretations of actual surge measurements. 
While dynamic magnetization measurements are to be preferred, 
nevertheless static flux data have proven very useful. Determination 
using such a magnetic circuit approach is illustrated in Fig. 13. The 
two curves are actual demagnetization curves (taken statically) for 
two typical telephone relay structures: a standard wire-spring relay, 
and a flat-spring miniature relay. For the particular operated conditions 
shown, the resulting values of (P are seen to be 1.8 X 10~ 8 Wb per 
ampere-turn in both cases. While these curves were obtained by measure- 
ments with a search coil, they could also have been estimated by well- 
known techniques of magnetic circuit analysis. 

The effective inductance may be inferred based on measurements, 
in a known circuit, of the transient behavior of such properties as 
frequency, critical damping, peak voltage, initial slope, etc. A few 
of these are noted below. 

3.1.2.1 Frequency Method. When circuit conditions are appropriate 
for an oscillatory response, as when G 2 <P/C and R t are approximately 
equal (discussed in Section 2.2), the frequency may be related to the 
circuit parameters based on relations given in Table II. Thus, 




100 200 300 

AMPERE TURNS, NI 



Pig. 13 — Demagnetization curves of typical relays. Upper: AF relay. Lower: MA 
relay. 
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/-" 



2tt 2irG, <P V/3i(l + Cj/G'x) 

1 

~ 2-jtN, VC(P(1 + 6' 2 /G,) 



or 



* ~ ktmcfQ + 0,/Qd' (28) 

Measured values of /, together with known values of iV, , G'i , and G 2 
are substituted in equation (28) to find the effective value of (P. 

3.1.2.2 Method Based on Initial Slope of V c . From Table V, note 
that when t = 0, V r = 0, for a load shunted by an RC network, when 

CR; - L,(l + G 2 /G w ) 2 

or 

65 = JV?(1 + (VG^ (29) 

Ordinarily, G 2 /G w is less than 0.01 and is neglected. Thus, in a circuit 
like Fig. 2a, one may vary R n and C in various ways until oscilloscopic 
measurement of the contact voltage surge gives a zero initial slope. 
Because of its convenience, this method has been used extensively 
in the present study. It also provides some picture of how the inductance 
varies with the shunt capacitance which is a rough measure of the 
speed of the transient. The method is illustrated in Fig. 14. Figure 14a 
shows a family of surges for different values of R n , with C fixed, and 
Fig. 14b shows a group of curves for which C and R n have been adjusted 
to give V c = at I = 0. The above procedure finds the combination 
of R n and C for zero initial slope, and substitutes these values, together 
with Ni , into equation (29). The results for a variety of representative 
armature-type relay structures are summarized in Fig. 15. It is seen 
that varying the direct shunt capacitance C, has a much larger in- 
fluence on the variation of (P, than does varying C in the RC network, 
presumably because the pulse transient is so slowed that it encompasses 
a larger range of the demagnetization curve. The spread in the data 
for various structures is thought to be explained by the range of con- 
ditions for magnetic structure, stop disc height, operated ampere-turns, 
etc. Figure 16 gives a summary of the effective permeance for a wide 
range of sealed magnetic contact relays; the value of (P is seen to be 
proportional to the number of contact units, as expected from magnetic 
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._* ^ |-a; 



;b) 



ILLUSTRATES VARYING SHAPE 
OF V c VS. t CURVE AS R n 
IS VARIED (C= 10" 7 f ). 

READING UPWARDS AT t = 0, 
VALUES OF R n ARE : 0, 500, 
1000, 2000, 5000, 10,000, 

100,000 ohms. 



ILLUSTRATES CONDITIONS 
FOR V c = AT t = 0. 

READING UPWARDS AT t=0, 
PAIRED VALUES OF R n AND 
C(fif) RESPECTIVELY, ARE: 
0, 0.1; 2700, 0.1; 2000, 
0.2; 1700, 0.3; 1500, 0.4; 
1230, 0.6; 1100, 0.8. 



Fig. 14 — Illustration of method of measuring effective inductance, based on V c = 
at t = 0. [AJ12 Relay; Ni = 5150; R w (hot) = 877 ohms.] 



circuit analysis. Each contact adds a permeance of about 0.12 X 10" 8 
Wb/NI per unit. Departures from complete "proportionality are at- 
tributed to different amounts of magnetic material in the return paths 
for various of these structures. 

The above values are found to check well with values for (P which 
are found as a by-product of other tests, as will be seen below. 

One may estimate the initial effective stored inductive energy, 
which will be released into the system upon disconnection, by means 
of the relation 

E = |<P(1 + G 2 /G W )(NI„) 2 (30) 

where NI e is the value of ampere-turns at the point of release. 

3.1.3 Method of Determining Effective Eddy Current Conductance, 
G e , oj Relay Core 

The properties of a coupled secondary winding may be calculated 
directly from the relation G 2 = N\/R 2 or from standard conductance 
formulas for copper sleeves. 4 When sleeves are used, their conductance 
values ordinarily outweigh the core's eddy current conductance to 
the point that it can be ignored. However, especially for unprotected 
applications, an estimate of the core's effective eddy current value is 
needed; and at the present time only experimental procedures seem 
suitable. Two such methods will now be discussed, one based on readings 
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of current or surge voltage at / = 0, and the other based on measure- 
ments of critical damping involving the use of a supplementary sec- 
ondary winding. 

3.1.3.1 Method Based on Value oj i t or V, at ' t = 0. An examination 
of the circuit variables indicates several equations where the term 
G 2 /Gi appears. For example, in Table V, the initial values of i x or V c 
involve the term 



K = 



1 ± G 2 /G w 
1 + G 2 /G, 



i\ 


at t = 


+ 


v c 


L 

at t 


= + 



Thus it is conceptually possible to measure G 2 by observing the value 
of K, experimentally, as by oscilloscope, for particular values of R n 
and R w . Rearrangement of the above expression then gives 



G 2 /G w = 



1 - KR, 
However, for small values of G 2 , the value of K will evidently be close 




6 10" 7 



C(f) 



Fig. 15 — Values of (P determined by V e = Oat I = method (conventional relays). 
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Fig. 16 — Values of (P determined by V c = at I = method (magnetic reed 
relays). 



to unity, requiring experimental accuracies that seem quite difficult; 
therefore this method seems of doubtful value for measuring small 
values of G 2 , though probably useful when G 2 /G w exceeds about 0.1 
as for copper sleeves, etc. 

3.1.3.2 Two-winding Method. Instead, the method using two windings 
as already derived in Section 2.2.1 seems more useful. The procedure is 
to vary the shunt capacitance across the primary winding (with R„ = 0) 
and observe the value of series resistance across the secondary (or 
vice versa) for which critical damping is recognized on an oscilloscope. 
These steps are repeated over a range of values for C in the region 
of 10" 10 to 10" 9 F. After correcting C for the distributed capacitance 
of the primary winding, the resulting values of 



G. = 



m 



R 2 + R' 2 



are plotted against C*. In accordance with equation (9), the data tend 
to fall on a straight line whose zero intercept represents G e . The effective 
permeance is then determined from the slope as 



9 - (&• <32) 

Figure 17 illustrates this procedure for a typical wire spring relay. 
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4N. 2 


2000 


(SLOPE) 2 




= 2.6x to -8 Wb/NI 




SLOPE = 8.1 X10 7 — ^/^ 


1000 


J3 





/^^- G e = 220 mho 

i i 



1 X10" 



Fig. 17 — Experimental determination of eddy current conductance, G t 



The data are from the small values for C and G 2 , in Table VIII, which 
will be discussed in connection with more extended experiments on 
critical damping (Section 3.2). In the present experiment, the inferred 
value of G„ , which is the negative intercept on the G 2 axis of Fig. 17, 
is found to be 220 mho, and the slope is 0.81 X 10 8 mho / _i , which 
from equation (32) gives <P = 2.6 X 10~ 8 Wb/NI. 

3.2 Critical Damping Measurements 

This section will describe measurements to verify the occurrence 
of critically damped relationships derived in Section 2.2, and to confirm 
the analysis for windings which are imperfectly coupled. 

3.2.1 Critical-damping Boundaries 

Determination of the effective eddy current term G t described in 
the previous section is a by-product of an experiment to verify the 
analysis of Section 2.2 which derived relations for the boundaries 
between oscillatory and damped conditions. It was particularly in- 
teresting to verify the case for core losses dominant, as this relation 
had not previously been recognized. The procedure, as previously 
noted, was to interrupt the primary winding of a two-winding test 
coil with a mercury contact (to minimize arcing) and to observe the 
character of the interrupted surge as the values of G 2 and of the ca- 
pacitance C across the primary were varied. (G 2 is varied by varying 
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the value of R 2 across its terminals). The value of the variable term 
for which the pulse changes from damped to oscillatory is fairly easy 
to recognize, using an oscilloscope. A set of such measurements is 
given in Table VIII together with calculations for plotting the data 
to check the validity of equation (5). The procedure was to (i) determine 
the corrected value of G 2 from equation (9), using the steps noted in 
the previous section and Fig. 17, (w) tabulate G 2 /G w , and (Hi) estimate 
the corresponding value of /3„ . 

A special consideration applies to this calculation of j8„ . The analysis 
in Section II has shown that the effective inductance increases with 
G 2 /G w (see Table I) as follows: 

V = (P'(l + G 2 /G w ). 

Here a prime has been added to the permeance term to indicate that 
(?' represents the actual state of magnetization of the coil's magnetic 
circuit. The value of (P in the formula ft, = CRI/&N] must be related 
to this experimental value as follows 

<P = <P'/(1 + G 2 /G w ). 

Thus, for the present measurements, a correction has been applied 
so that py, can be evaluated using the experimental value, (?'. The 
corrected value is 

p w = (CRt/V'N\)(\ + G 2 /G w ). (33) 

The values of W as estimated in this way, and as read from the curve, 
of Fig. 4, are listed in adjacent columns of Table VIII, and the experi- 
mental points are plotted as circles in Fig. 4. Agreement between 
experiment and analysis is seen to be good, especially considering 
that a constant value of (P' has been assumed to apply over the whole 
range of the experiment. 

The family of curves in Fig. 5 provides a convenient method for 
analyzing the data for a particular structure. In this case one plots 
the value of G 2 /G w vs the corrected shunt capacitance at critical 
damping. Then, based on the known values for JV, and R w for the 
given structure, one infers the effective value for (P' from the numerical 
value for the appropriate curve, using the correction technique of 
the preceding paragraph. Data for several different structures with 
secondary windings are plotted, and also for several relays with copper 
sleeves. The copper sleeve cases are identified by a dot (measured 
value) surrounded by a rectangle representing the possible spread 
in the various parameters for a sample relay drawn from stock. The 
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copper sleeve relays represent the four standard sizes used in slow 
release relays for telephone switching. Agreement is again seen to 
be good. 

3.2.2 Oscillations in Im perfectly Coupled Circuits 

Measurements using secondary windings give somewhat different 
results depending on the form of winding used. When the primary 
and secondary are perfectly coupled, as with parallel-wound coils, 
the boundary between oscillatory and damped regions accords with 
the analysis given above. However, when primary and secondary 
are not perfectly coupled, as in the case of concentric windings, or 
of copper sleeves, a new and at first unexpected result is noted. In 
this case, there is a new form of oscillation which becomes noticeable 
as one operates near or beyond the region of critical damping — this 
is the discharge of the uncoupled portion of the inductance which 
is not attenuated by the core-loss term. The analytical basis was 
discussed in Section 2.2.2 and some data on numerical values will now 
be given. 

The presence of shorted secondary windings which are not perfectly 
coupled to the primary winding was seen to lead to two superimposed 
transient patterns: (t) that due to the coupled portion of the winding 
which experiences critical damping below values of determined by 
the coil losses, and (it) that due to the uncoupled portion of the winding. 
For this group of windings, operating in the overdamped region, one 
observes oscillations in the circuit which are produced by the parallel 
inductive path of the uncoupled part of the main winding. Such oscil- 
lations are particularly noticeable with relays using copper sleeves. 
In fact the experimental observation of such oscillations is often helpful 
in the experimental determination of the critically-damped condition 
for the coupled portion of the winding. The following results are given 
as guides for dealing with these effects, which often need to be con- 
sidered in the course of relay application studies. They are taken from 
extensive measurements on copper sleeve relays. 

Among the general results were (?) the consistent observation that, 
when the coupled portion of the main winding was operating in the 
core-loss overdamped region, the amplitude and frequency of the 
oscillations for the uncoupled portion of the coil were independent 
of the value of G 2 , and (ii) the frequency was proportional to C , 
as predicted by equation (11). In addition, the critical damping bound- 
ary for four sizes of sleeve (covering a range of G 2 between about 
40,000 and 200,000 mho) was measured. The data were then analyzed 
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to give the corrected value of effective core permeance for different 
values of Q*/G\ , in the same manner discussed in Section 3.2.1. This 
value of permeance (P' is determined by using the value of /3,„ for a 
given value of G 2 /G w from the curves, Fig. 4, and substituting in 
equation (33), as is shown in the left-hand side of Table X. The spread 
in the indicated value of <?' seems reasonable for a group of randomly 
selected stock relays. 

The right-hand portion of the table shows measured values of fre- 
quency when the primary was shunted by a capacitor, C = 10~ 8 F. 
Now, calling the effective value for this uncoupled inductance L" = 
(P"N'\ , all the test samples are seen to show values for <P" of about 
0.6 X 10" 8 Wb/JVJ. This value can be used to estimate /3, V p , /, etc., 
for the uncoupled portion of the winding, for which G 2 = 0. 

These data indicate that the uncoupled portion of the main winding 
can contain an appreciable energy which is externally releasable, 
even though the core-loss term due to the copper sleeve is very large. 
For example, for the relays studied in Table X, the initial releasable 
energy is about E = 0.6 X 1Q~*(NI ) 2 joule. 

The above concepts explain the general transient behavior of copper 
sleeve relays. The numerical values which are given for (P only apply 
to the particular structural group of relays studied here (the slow- 
releasing class of the wire spring family). 

3.3 Time-to-Peak, t', and Peak Voltage, V p 

Peak voltage measurements, i.e., measurement of the time-to-peak 
and the value of the peak voltage, offer further opportunities to check 
the analysis. The following discussion examines the measurements 
for both peak time and peak voltage on a single experimental unit, 
and also reconsiders a large series of peak voltage measurements which 
were made about 1953 based on then-current empirical practices. 

3.3.1 Measurements of t' and V p on a Sample Wire Spring Relay 

The values of t' and V p were measured on a test relay over a range 
of shunt capacitances, C, representing a large spread of the variable /3, r . 
Based on the known constants of the structure, the formula for t'/R w C, 
equation (12) as plotted on Fig. 7, and discussed in Section 2.3, was 
used to infer the associated value of (3„. . This value of 0„, was then used 
to predict the expected value of V p /V n , from Fig. 7. These results are 
summarized in Table XI, together with the measured values of V p /V„ . 
The measured and predicted values for VJV are seen to agree very 
closely, and are taken to support the analysis. 
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The final column in Table XI indicates the effective value of (P for 
this same series of measurements. It is seen to have a fairly constant 
value for values of fi w greater than about 10~ 7 , and to gradually in- 
crease above this point. In line with earlier discussions, this increasing 
trend accords with the view that slower surges involve larger changes 
in the demagnetization curve and correspondingly larger values of (P. 
In keeping with this picture, it appears that the values of (P as deter- 
mined at the initial stages of the pulse (from V B — at t = in Section 
3.1.2) should be increased by a factor of 2-3 when working with events 
at the peak of the pulse, where the flux curve has had more time to 
change. 

3.3.2 Use oj 1952 Data to Further Confirm Peak Voltage Relation 

During a two-year period starting in about 1952, an extensive study 
was originated by M. M. Atalla, and later extended by many others, 
into the values for peak voltages when contacts open relay loads. 
The data were found to fit a formula relating VJV, to a variable 
similar to the W term used in the present study, modified by an ex- 
ponential term that was determined empirically. These very extensive 
data are available for comparison with the present analysis; their 
sources are unpublished Bell Laboratories memoranda. This earlier 
study used the following notation for the term designated /3„ in the 
present work: 

where (R„ represented the reluctance of the relay's magnetic circuit, 
and is the reciprocal of (P, the permeance as used in this report. This 
reluctance value, technically described as the "closed-gap reluctance," 
is a design constant which is tabulated for all telephone type relays. 
In the light of the present analysis, this value of reluctance, which 
represents the steepest slope of the magnetization curve, should be 
replaced by the uppermost slope of the demagnetization curve, repre- 
senting a much higher reluctance (or lower permeance). Thus, all the 
earlier data can be used as a check on the present analysis if a relation- 
ship between these two figures is known. Data are not available for 
all the specific cases involved; however, it was determined from a 
review of magnetic data for many of these structures that a correction 
from the (R values mentioned above to values for (R which are related 
to the demagnetization curve could be approximated by increasing 
(R„ by the factor 7.5. When this is done, when units are changed from 
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cgs to SI, and the new values of /3„, are calculated, and (R is converted 
to (P (its reciprocal), a large number of data points for different struc- 
tures becomes available for use here. The revised data points, taken 
from the unpublished memoranda, are plotted for various structures 
on Fig. 7. They are seen to show very good agreement between analysis 
and experiment for a wide range of telephone relays. The agreement 
could be further improved at the longer transient times (larger values 
of C or j3,„) if an added correction for increasing differential permeance 
were made, along the lines indicated in Table XI (column showing 
"Inferred (P.") 

Although the core-loss terms for these particular structures have 
not been measured by the methods outlined in this paper, the data 
points as plotted for the various structures in Fig. 7 are consistent 
with the expectation that large cores, or low-resistivity cores, should 
have larger values of G a /G w than small or high resistivity cores. These 
results therefore offer additional confirmation for the analysis of Section 
II. 



SECTION IV. DISCUSSION 

Based on the close agreement between analysis and experiment, 
the formulas given here appear well suited for use in analyzing the 
various transient effects associated with a contact opening an inductive 
load. Because of satisfactory prediction of pulse times, peak voltages, 
frequencies, and regions of damping or oscillation, it also seems reason- 
able to extend the use of these relations to the approximate estimation 
of energy dissipated in the contact arc, according to the analysis of 
Section 2.4. 

Returning to the list of various problems discussed in Section I, 
the present study should therefore find application in the following 
areas. 

(i) Prediction of surges from relay windings as a function of time, 
including breakdown values, times, oscillations, etc. The results 
apply equally to loads controlled by transistors, vacuum tubes, 
or relay contacts. While the examples are mainly given for 
the "unprotected" load, the basic formulas are available for 
the study of other cases. 
(it) Extension of the above work to predict the energy which will 
be internally dissipated or be made available via the capacitor 
to do damage in arcing, as in contact erosion. Extension of the 
present work is still needed to more fully evaluate the energy 
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replenished into the system from the battery during arcing. 
It is expected that continuing work" will improve understanding 
of the present empirical technique which relates contact erosion 
to the sum of an energy term and a current term. 
(in) Reasonable explanations now are available for earlier uncertain 
points. 

(a) Relays with higher stop discs (i.e., larger operated air gaps) 
will ordinarily be less close to saturation. Therefore, their 
differential permeance will be larger than otherwise, which 
predicts a larger releasable stored energy. Although this 
configuration does not represent the condition for maximum 
total stored energy, (which must be at the smallest air gap), 
it does represent the condition for large initial release of 
magnetic energy, i.e., the value of (P is larger. 

(b) The oscillations which occur with copper sleeve relays are 
found to result from the uncoupled portion of the main 
winding. Relations are given for predicting their values. 

(iv) Methods for predicting various behavior patterns are demon- 
strated. 

(a) An example of a "B" type transient energy calculation is 
given. 

(b) An extension of (a) above can be used to predict the final 
peak surge, following the last breakdown of a "B" type 
transient. Particularly with coils having a large number of 
turns, and a slowly rising surge, the number of breakdowns 
in a "B" type transient is relatively few. In such a case, 
a considerable residual inductive energy may remain 
following the last arc breakdown. The associated final surge 
tends toward higher than ordinary peak surges which need 
to be evaluated in order to insure adequate protection 
against breakdowns in the associated wiring. 

(c) An analytical basis is available for determining the proper 
test conditions for coils being production-tested to with- 
stand their own self-induced surge voltages. 

(d) Many past evaluations of contact wear, and comparisons 
of various alloys and environments, are hard to compare 
with metals being studied today because the test conditions 
(magnetic structure of the load, circuit voltage, etc.) were 
so different from today's conditions. The relationships in 
this report are expected to permit the re-evaluation of some 
of this older data in a form that will permit comparisons on 
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a more universal basis. In one such test (1923), more than 
250 metals and alloys in various gaseous environments 
were tested and compared. The review of this information 
in the context of today's applications problems should 
help to guide current design activities, as for example, 
those on sealed contacts. 
(v) Future design for relay miniaturization can be better evaluated 
based on fuller understanding of the influence of the inductive 
load on the wear (and hence the necessary volume, travel 
allowance, and required work) of the contacts which control it. 
For example, the basic element determining contact wear is the 
effective releasable stored magnetic energy. The principal factor 
in this quantity is the initially stored effective energy term 
%(?(NI )' 2 , and often an additional term proportional to current 
/„ due to energy replenishment into the system by the battery, 
during arcing. Thus, contact erosion is related to the square 
of the operated ampere-turns and the differential permeance 
of the magnetic circuit's demagnetization curve, together with 
a term of approximate value 0.1 /„ . The most favorable con- 
dition is therefore that where the desired work output and 
operating time for the relay lead to a low value for NI (and 
also of current). The low value of NI is generally associated 
with a high-permeance magnetic circuit at the unoperated or 
critical gap condition. A high permeance at the closed-gap 
releasing current, however, is not compatible with low re- 
leasable energy. The present study suggests that contact life 
can be further improved if this sensitive unoperated structure 
were designed so that it is highly saturated (i.e., low (?) once 
it is in the operated position. 

This releasable energy may be further reduced through the 
conscious use of a secondary winding whose value of G 2 is large 
enough to decrease the energy available for arcing by several 
fold while not being large enough to seriously influence release 
time. Typically, this means that one might profitably use a 
small closed secondary winding or copper sleeve having a value 
of G 2 equal to about 1000 to 5000 mho. If desired, the effect on 
operate time could be minimized by completing the secondary 
winding circuit through a diode which is poled to make the 
secondary effective only on release. 

If, in this way, the load magnet can be made to cause less 
contact wear, the smaller contact wear allowance can result in 
a lowered mechanical work requirement. This in turn leads to a 
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smaller magnet structure, and so on. The concept should thus 
prove useful in design steps aiming to optimize a miniature 
relay structure, including the sealed magnetic reed types with 
their thin electroplated contact regions. 

V. SUMMARY 

This study has presented a group of analytical relations which 
describe the electrical transients when a contact opens an inductive 
load. In the past, relations based on a simple LRC analysis were found 
to be quite inaccurate. The present method expands the earlier approach 
to account for core losses, and is found to reduce to LRC-type relation- 
ships provided the lumped-constant terms have effective values as 
defined in Section II, and provided the value used for inductance 
recognizes the actual differential flux-ampere-turn characteristic for 
the particular load's magnetic circuit. 

The analysis then leads to a group of useful relations for predicting 
surge times, peak voltages, or energies either consumed internally 
or available for external damage as in contact arcing. The analysis 
also predicts the conditions under which the surges are either oscil- 
latory or damped. Experiments have confirmed these results, including 
the presence of two critically damped regions whenever the load pos- 
sesses a finite secondary loss term. They have also explained the presence 
of an "anomalous" oscillatory condition for two-winding structures 
resulting from the uncoupled portion of the primary inductance. 

Illustrations of how these results may be applied in circuit and relay 
design are given, as for example: the prediction of peak voltages, or 
of the chain of breakdowns in a "B" type transient; the estimation 
of energy available to produce contact erosion; and the determination of 
the oscillation frequencies in a given circuit. Experimental techniques 
were also described for determining the effective values of the various 
circuit constants, and typical values for a wide range of structures 
are given. The implications in relay design, as for miniaturization, 
are also discussed. 

APPENDIX A 

Derivation of Basic Circuit Equations for Current 

The Disconnection of a Relay Winding-General Case (Assuming Ideal 
Break Contact, an RC Shunt Network, and a Secondary Winding) 
Upon disconnection, transient currents will flow in the two meshes 

as shown in Fig. 2a. Using the notation of that figure, and of Table I, 
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the equations are 



Nl Tt + i,Rl + 



/* 



dt 



C 



= 



= 



(34) 



These are interrelated through the magnetic circuit, for which we 
assume * = 3>, + 4> 2 , where <f>! and <f> 2 are proportional to the primary 
and secondary magnetizing forces, N,i, and N 2 i 2 , respectively. Now 
designating (P as the incremental permeance of the relay's magnetic 
circuit, i.e., (P = dip/dNi, at the point in question, we write 



d± JdNjh , dN,i, 
dt == \ dt + dt 



(35) 



Making this substitution, the loop equations become 



N JdNA + dNA ] + iiRi , 



dt 



dt 



h 



</t 



C 



v^ + ^) + - 



= o 



= 



(36) 



We now adopt LaPlace Transform notation, and combine with initial 
conditions, which are that at t = 0, the charge on C is ( — V ), and 
the core flux is $„ = (?NJ n . Equation (36) may then be rewritten: 



NtfaiNj, + N 2 i2) 4- hB, + -^ = NtfNJo + ^ 



NiVsiN,!, + iV 2 z 2 ) + uR> 



= iV 2 (PiV,Z 



(37) 



Using the second of these equations to find i 2 , substituting into the 
first, and collecting terms: 



h = L 



NMl + (! 2 /a, r )s + R, 



Nfra + G 2 /G,)s 2 + [r, + ^)s + ^ 



(38) 



This is recognized as the expression for a simple R, L, C circuit, as 
shown in Fig. 2b, for which the inductance and resistance have effective 
values: 
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U = L,(l + G 9 /G x ) = Nl<P(l + G 2 /G x ) 
R' = R x + G 2 (P/C. 
Using these values, equation (38) becomes 

(s + a) 



ix = KI 



(s + b) 2 + 



(39) 



where values of a, b, and to are as given in Table II. Equation (39) is 
immediately solved by reference to tables of LaPlace Transforms. 
There are three cases: 

I. Oscillatory: io 2 > b 2 

II. Critically Damped: u 2 n = b 2 

III. Overdamped: of < b 2 

The solutions are found to be: 

Case I. 

KI e~ bl sin (wt + <p) 



ii = 



where <p = tan 



sin <p 

-1 w__ 

a — b 



Case II . 



Case III. 



t, = KI„e~ bl [l + (a - b)t] 
KI e- b ' sinh (tat + <p) 



i\ = 



where <p = tanh" 



sinh <p 

CO 

a — b 



(40) 



(See Table III for more detailed summary of values for y and related 
terms.) 



APPENDIX B 



Evaluation of Maximum or Minimum Values of q x , i x , and V e 

The procedure for these determinations is to differentiate the cor- 
responding expression for q x , i t or V c given in Table IV, with respect 
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to t, and set the result equal to zero. This result, which defines the 
times at which the slope of the transient = 0, then determines either 
a maximum or a minimum condition, as will be shown. The first maxi- 
mum or minimum will be designated t' , i.e., "time-to-peak". The 
value for t' is then substituted into the original expression to give the 
"peak" value for q x , i t , or V c (using the subscript p to designate it). 
While the analysis shows the steps taken for the general case of a shunt 
RC network, the resulting equations are often fairly cumbersome. 
Thus, after indicating the direction of these general results, they are 
fully expressed for the case of shunt capacitance alone (R„ = 0) which 
is the major subject of the present study. The steps are briefly outlined 
below: 
Case I — Oscillatory (i.e., u 2 > b 2 ) 

Charge 

From Table IV, the time derivative of charge is 

% = ^- | le-"' sin M + + *)] 

at u„ sin <p at 

A ' "*> cos (at + 6 + <p) - b sin (ut + 6 + <p)\ - 0. 



oj sin <p 

The derivative thus has a zero value when t' = °o or, based on phase 
angle relations indicated in Table III, when 

tan (at' + 6 + <p) = tan 6, 

which is only satisfied at values of /': 

l " ' o, J o, ' GtC - (41) 

(max) (min) etc. 

Substituting equation (41) into the original expression for q Y : 

?,„ = +^~ sin 6 exp (-^Fr ) ■ (42) 

* v <j) sm<p r \ tan 5 / 

For the case when R n = 0, from the relations in Tables II and III, 
this reduces to 

ft.-CF*»-p(-^f)- (43) 

Current 

From Table IV, the time derivative of current is 



TRANSIENT EFFECTS 



3057 



di x Ki d . _(,, . / , , \i 

» = r. le sin (co/ + <p)\ 

dt sm <pdt L v J 

= - — -e bi [u cos (co£ + <p) — b sin (coi + <p)]. 
sin ^> 

The derivative thus has a zero value when t' = °o or, based on phase 
angle relations of Table III, when 

tan (ul' + ip) = tan 0, 

which is only satisfied for values of t'\ 

_ 7T + - ^ 2ir + - ip 



t\ = 



CO 



mm 



max' 



, etc. 



, etc. 



(44) 



Substituting equation (44) in the original expression for U : 
Ki„ sin 



sin tp 



Or + - ip) 
tan 



For the case when R n = 0, from the relations in Fig. 4, 

Or + - ip) 



Vl 4- Go/G w 



- exp I — 



tan 



(45) 



(46) 



Contact Voltage 

From Table IV, the time derivative of contact voltage is 



dV e Ki„ d 



dt sin ip dt 



— e b, \R v sin (ut + ip) ^ sin (at + <p - 



C 



= -^ e~ bt \ ^sin (wd + fp) - wje,, sin (co* + <p - 0) 
Bin <£> L^ 

The derivative has a zero value when /' = °° or, based on phase angle 
relations of Table III, when 

tan 



tan fat' + <p) = 



which is satisfied for values of t'\ 



1 



1 



Cu n R n cos 



CO 



tan" 



tan 



1 



CwJ2, cos 



(47) 



* (usually, however, there will be an initial maximum preceding min*, at (x/2 — 
y), when >f < 7r/2). 
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The general expression for peak voltage is then given by substituting 
equation (47) into the original voltage expression. The resulting rather 
complicated formula may be used to estimate peak voltages for the 
general case which includes an RC protection network, a procedure 
which is beyond the scope of the present report. For the case when 
R n = 0, substituting values of the various terms from Tables II, III, 
and IV: 



/' = --- -- . ^LZL* . etc. 



IT — <p 



max , min , etc. 

(48) 



V cp = V.[l + «? exp ( -* ~ s 



tan 9 

Case II — Critically Damped (i.e., u>l = b 2 ) 

For the critically damped case, as for Case I, the procedure is to 
find the time, t', at which zero slope of the time derivative of each 
variable occurs, and the corresponding "peak" value, for the general 
case which includes the RC network. Following this, the results are 
given for the specific case for R n = 0, with individual consideration 
of the two critically-damped conditions: core losses dominant, and 
coil losses dominant. 

Charge 

From Table IV, the time derivative of charge is 

= -Ki e~ bl [(b - a)t - 1]. 
The derivative has a zero value when t' = »or when 

t' Q = 1 - } ~- (49) 

b — a 

Substituting equation (49) into the original expression for q: 

CV \KR' J / b \ ,_ m 

or 

= 0. 
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Current 

From Table IV, the time derivative of current is 

= -^>-"[6(6 - a) - 26 + a]. 
The derivative has a zero value when t' = °o , or when 

'' - ^T=^)' (51) 

Substituting equation (51) into the original expression for ?', : 

„. (6 — a) / 26 — q \ . . 

„., = -Ku — — exp ^-y^J 0-2) 

or 

= 0. 

Contact Voltage 

From Table IV, the time derivative of the contact voltage is 
dV v d[ _„/ KB. , (b KR\\ 

~dT = - v °itV i 1 - ~r~ + (a - 6) la - urn 

Ki 
= ^e""[l - CR„(2b - a) + (a - &)(1 - WR„)t]. 

The derivative has a zero value when t' v = <x> , or when 

_ 1 - Cfl„(26 - a) 
r (6 - q)(l - />Cfl„) V } 

and its value may be found by substituting equation (53) into the 
original voltage expression. 

Cases Where R„ = 

A more detailed treatment for the case where R„ = will now be 
given, to indicate behavior under the two conditions for critical damping: 
(0 Core losses dominant: i.e., 

( -~ > R, r , or ^ < /3, , or b > a (0, > 0) 

(?'/) Coil losses dominant: i.e., 

R„ > <J ~ , or /3, r < ^ , or a > b (/3„. ^ 4) 
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Then, the values of minimum or maximum circuit variable, and of 
corresponding times, are summarized below: 



(i) Core Losses Dominant 
Charge 

1 



ti = 



b — a 
= <7,„(P/3i(l + A) 2 



q u> = CV ($ 



~i„- (1 



Current 



t' = 



26 - a 



b(b - a) 



Contact Voltage 

i 



• 2)03* + 1) 



b — a 



= r;„.(p/3?,.(i + oi) 2 
v cp = F [i + ^ e - (1 ^"; 



(it) Coil Losses Dominant 



ti = - 



a — b 



As a > b, a negative, i.e., 
impossible, value of /' is in- 
dicated. Thus, t' q = qo applies, 
and there is no peak value of q, 
which gradually decays to 0. 



t' = - 



2b - a 

b(a - b) 



As just above, this expression 
is always negative, and t' t = °° 
applies, indicating no mini- 
mum, but a gradual decay 
toO. 



ti = - 



a — b 



As above, this expression is 
always negative, and t' v = » 
applies, indicating no maxi- 
mum, but a gradual rise in 
contact voltage from to V„ . 



Case III — Overdamped (i.e., 6" > to 2 J 

By following steps similar to those for Case I, except noting that 
the equations are in hyperbolic rather than trigonometric form, the 
conditions for maxima or minima are found to be as follows: 

Charge 

The condition for a maximum is that t' = °o , or 
tanh (uf + 6 + ip) = tanh 0. 
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The latter case is seen to be satisfied when 

< = -I- (54) 

Now, recalling that <p = tanh -1 co/(a — b), it is seen that a positive 
(i.e., realistic) value of t' only occurs when b > a, which is the case 
where core losses are dominant. Therefore when coil losses are dominant, 
the charge gradually decays to zero. 

The resulting peak value of charge, when b > a, is 

Ki n sinh 6 ( <p \ ,__. 

q >» = co.sinh* eXP \ _ tlmhW (55) 

where cp = tanh -1 o>/(6 — a). 

For the case where R n = 0, this reduces to 

9 ,„ = CK„ft'exp(- tl ^). (56) 

Current 

The condition for a minimum is 

tanh (cot' + <p) = tanh 5 
which is seen to be satisfied when /' = oo , or when 

CO 

As <p is always greater than 6 (see Table III), /' can only have a positive 
value when b > a, for which case 

6 + tanh - ' " , 



t' = -J^-*. (57) 

CO 

There is no minimum for the case that a > b, except as the current 
gradually decays to zero. 

The resulting peak value of current, when b > a, is 

KI n sinh 6 ( 6 - y \ ,„. 

For the case where R„ = 0, this reduces to 

{ " - Vrta7<r„ exp ( - 4ffi#) ' (59) 
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Contact Voltage 

The condition for a maximum is that i' = °° , or that 

tanh K + <p) = ^~ (60) 

1 

Ca) R n cosh 

This expression may be used to give the peak voltage, when b > a, 
for the general case, as was indicated above for the oscillatory case. 
For the case R n = 0, the primary subject of the present paper, this 
becomes 

tanh (cot' + e) = 
or 

t' - ■=*• (61) 

CO 

As in the case of equation (54), this indicates a peak value of voltage 
only for the case where core loss dominates, for which 



V = V 



1 + p? exp - 



(62) 



tanh 6 1 
When a > b, the voltage gradually rises to the asymptotic value V c = V„. 

Summary 
The foregoing expressions are summarized in Table VIII. 
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